Monte Carlo Study of a 17(1) x U(l) Loop Model with Modular Invariance 
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We study a U(l) x 17(1) system in (2+l)-dimensions with long-range interactions and mutual statistics. The 
model has the same form after the application of operations from the modular group, a property which we call 
modular invariance. Using the modular invariance of the model, we propose a possible phase diagram. We 
obtain a sign-free reformulation of the model and study it in Monte Carlo. This study confirms our proposed 
phase diagram. We use the modular invariance to analytically determine the current-current correlation functions 
and conductivities in all the phases in the diagram, as well as at special "fixed" points which are unchanged 
by an operation from the modular group. We numerically determine the order of the phase transitions, and 
find segments of second-order transitions. For the statistical interaction parameter 9 = n, these second-order 
transitions are evidence of a critical loop phase obtained when both loops are trying to condense simulataneously. 
We also measure the critical exponents of the second-order transitions. 



I. INTRODUCTION 

Models with statistical interactions can be used to describe 
a variety of interesting systems. In particular, quasiparticles in 
the Fractional Quantum Hall effect, as well as other fraction- 
alized phases of spins and bosons, have such interactions^^ 
Some models with statistical interactions contain a symmetry 
under the action of the modular group. This can simplify an- 
alytic study of these models. Several different such systems 
have been studied in the literature^— In this work we define a 
model with this symmetry, which we call modular invariance, 
and study its properties both numerically and analytically. 

In this work, we study a U ( 1 ) x U ( 1 ) model in (2+ 1 ) dimen- 
sions with mutual statistical interactions. After introducing 
the model, we will explain what we mean when we say that 
it has modular invariance. A general action for two species of 
U(l) particles with mutual statistical interactions is given by 
the following action: 
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J2i9Jx(-k) ■ a 2 (k). 




(1) 



Here J\ and J 2 represent conserved integer-valued currents 
residing on interpenetrating cubic lattices, and V • J\ = 
0, V • J 2 =0. For brevity, the above action is defined in terms 
of Fourier components, where v\(k) and v 2 {k) are Fourier 
transforms of the intra-species interactions for species Ji and 
J 2 respectively. In the partition sum, a given current con- 
figuration obtains a phase factor e lB or e~ %B for each cross- 
linking of the two loop systems, dependent on the relative 
orientation of the current loops. This is realized in the last 
term of Eq. (HJ, by including an auxiliary "gauge field" a 2 , 
whose flux encodes the J2 currents, J 2 = V x a 2 . As ex- 
plained in our previous worksriii the model is precisely de- 
fined with periodic boundary conditions if, for all directions 
fx, Ji M)to t = J2 r J it*( r ) = °> and similarly for J 2 . 

Of some relevance to our study is the work of Fradkin and 
Kivelson, - Though several of the mathematical results in their 
work are applicable to our model, the model itself is different. 
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FIG. 1: Phase diagram for the model in Eq. 0, with fixed 9 = 
2tt/3. The dashed line is the 'symmetric' line where the potentials 
are equal, vi = 1)2, which is also assumed everywhere in the phase 
diagram in Fig. [2] In phase "00" both J variables are gapped, while 
they condense in phase "0". In phases I and II only one species is 
gapped. In the lower left corner different composite variables are 
gapped; here the structure can be significantly different for different 
values of 9. 



In particular, Fradkin and Kivelson from the outset require a 
binding between different species, which is not present in our 
model and does not occur in our phase diagram. 

Figure Q] shows a schematic of the phase diagram of the 
model in Eq. £[), for 9 = In the remainder of this work, 
we use the following marginally long-ranged potential: 
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(2) 



where \f k \ 2 



^2 „(2 ™ 2 cos kfj,), and gi 2 are parameters de- 
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scribing the strength of the potentials. The v(k) ~ pry 

mentum space is equivalent to a 1/r 2 potential in real space. 
The main features of the phase diagram are controlled by the 
overall strength of the potentials, independent of the precise 
form (e.g. in a previous workii we found a similar phase di- 
agram for short-ranged potentials), while more detailed prop- 
erties do depend on the range of potentials. The dashed line 
in the figure is the 'symmetric' line, where vi = V2- 

When vi and V2 are large, both Ji and J2 particles are 
gapped, and only small loop excitations are possible in these 
variables. We call this phase "00". If we decrease V\ and v 2 
along the symmetric line, the reduction in potential allows the 
Ji and J2 particles to condense, in the sense to be defined 
later; we call this phase "0". The labels of these phases will 
be explained in Sec. [TTIJ Having both particles condensed at 
the same time incurs some penalty due to the statistical in- 
teraction, and phase exists only at intermediate couplings. 
As vi and v% are reduced still further, we can reach a phase 
where multiples of J\ and J2 (in particular, multiples of n if 
9 = 2ir/n) can condense. Roughly, such n-tuples of the cur- 
rent variables do not see a statistical interaction, the more pre- 
cise meaning of this will be explained below. Off the symmet- 
ric line, we can access phases I and II where only one species 
of loop is condensed, and the other is gapped. The phase di- 
agram is qualitatively similar for other values of 9, with two 
exceptions. First, for 9 = it phase is not present. Instead 
there is a phase transition on the symmetric line, which sep- 
arates phases I and II. An example of such a phase diagram 
can be seen in Ref. Q~H Second, for generic values of 9 more 
phases exist at small g. These can be seen in vertical cuts in 
Figs.|2]and[3] and will be explained in Sec. Hill 

The modular group is an infinite, non-abelian group, gen- 
erated by two operations: duality (denoted by S) and period- 
icity (denoted by T). We call our action "modular invariant" 
because it has the same qualitative form after the application 
of these operations. The periodicity operation corresponds to 
shifting the statistical angle 9 by an integer multiple of 2ir, and 
since the loop cross-linking number is an integer we can see 
that e~ s for the action in Eq. (fl]i is unaffected by such shifts. 
In what follows we will find it useful to define 77 = 4-, and 
the complex number 



in mo- the following "dual" action 
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z = rj + ig. 



(3) 



In this notation the action of such a shift by an integer n is: 



T n : z -> z + n. 



(4) 



Duality corresponds to performing a well-known duality 
transform^— on both species in the above action to obtain 



<Sdual[<3l,<32] = 77 [«l : dual|Ql(fc)| 2 + ?>2,dual | Q'2 (fc) | ' 
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+ }JOdua,iQi(-k) ■ aQ 2 (k), 



(5) 



^1/2, dual 
#dual 



(2ir) 2 v 2/1 (k) 
\f k \ 2 v 1 (k)v 2 (k) + 9 2 ' 

-{2w) 2 9 
\fk\ 2 v 1 (k)v2(k)+9 2 ' 



The Q variables are dual to the J variables and are also 
conserved integer-valued currents satisfying V • Qi = 0, 
V-Q 2 = 0. Under the exact duality Qi, to t = Q2,tot = 0; clq 2 
is an "auxiliary" field such that Q 2 = V x oq 2 . If we think 
of the J variables as boson number variables, the Q variables 
are vortices in the boson phase variables. 

Let us use the long-ranged potential in Eq. (0, then we can 
see that on the symmetric line the action (O has the same form 
as (Q]). The parameters in the original action transform under 
the duality in the following way: 



ffdual 
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?7dual — 

_|_ jyZ g + ?7 

In terms of the complex number z we have 

S : z -> -1/z. 



-v 



(6) 



(7) 



Transformations S and T generate the modular group of trans- 
formations of the upper half of the complex plane z. There- 
fore with this choice of potential we say that the system is 
modular invariant. What happens here is that the statistical 
interaction can also be viewed as a marginally long-ranged in- 
teraction, and the duality operation preserves the form of such 
interactions. We chose the potential in Eq. (f2]i for analytical 
convenience because its form is preserved under duality, but it 
also corresponds to a three-dimensional Coulomb interaction 
between charged particles constrained to two spatial dimen- 
sions, and so we can apply this model to realistic systems. 

The phase diagram of a modular invariant system can 
be determined entirely from the properties of the modular 
group^iii We will also use these modular transformations 
to characterize each phase of our model in terms of quasipar- 
ticles gapped in that phase. This will allow us to determine 
the behavior of current-current correlators and conductivities 
in each phase. 

Our numerical study also allows us to examine the critical 
properties of the system. All of the phase transitions in the 
modular invariant phase diagram can be mapped to each other 
under modular group operations. Therefore all such related 
phase transition points should have the same critical proper- 
ties. We have found some phase transitions which are second- 
order, with continuously varying critical exponents. This is an 
example of a novel type of phase transition. We have also 
studied special points in the r/,g plane where three phases 
meet, and found these to be first-order. 
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II. MODEL AND MONTE CARLO METHOD 

The action in Eq. (Q]) has a sign problem, which must be 
eliminated if we are to study it in Monte Carlo. In order to 
do this, we dualize only one of the two loop species. In this 
work, we dualize the J\ variables to get>ii 



(2tt) s 



v 2 (k) 



k L 
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\QiW\' 
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\fk\ 2 Vl (k) 



47TgQi(-fc)- J 2 (fc) 

\h\ 2 vi(k) 



This is a sign-free classical statistical mechanics problem in 
terms of closed lo ops Qi , J 2 and works for any v\,v 2 and 9 
(Note that in Refs. fl3|fl4 we used a different sign-free refor- 
mulation that only works for a specific short-ranged v\, v 2 ). 
In order to study the above action numerically, we write it in 
real space and use the potential from Eq. (fJJ: 



S[Qi,J 2 } = Jxv(r-o 



-Qi(r)-Qi(r'){9) 
9i 



+■ ^) Mr) ■ J 2 (r') + ^Qx(r) • J 2 (r') 



9i 



V(r 



1 2?r 



Jk-(r-r') 



(10) 



where Vol = L 3 is the volume of the system. In real space, 
J 2 fj, (r) is an integer- valued current on a link r, r + fi of a 
cubic lattice. The variables J\ are defined on a lattice dual 
to the lattice of the J 2 , but after the duality procedure the 
Qi are integer-valued current variables on links of the same 
cubic lattice as the J 2 . We perform our simulations using 
the directed geometric worm algorithm^ We attempt to pro- 
duce worms in both the Qi and J 2 variables, while satisfying 
Qi.tot = J 2 ,tot = 0. In this work, we monitor "internal en- 
ergy per site" e = 5*/ Vol, and compute "specific heat", defined 



C = ((e 2 ) (e) 2 ) x Vol. 



(11) 



In what follows, we will present all of our results in the 
Ji, J 2 language of Eq. (|T). To study the behavior of these 
variables we wish to monitor current-current correlations, de- 
fined as: 



C£(k) = (Ja»{k)J b „{-k)) 



(12) 



where a and b are the loop species and \i and v are direc- 
tions; J afl (k) = 77= Er J «m( 



= »*-r_ \Ye trivially have 

C~(fc) = C~,{— k). Because of the vanishing total current, 
we define the correlators at the smallest non-zero fc; e.g., for 
C%% we used k = (0, ^ , 0) and k = (0, 0, ^ ). For simplic- 
ity, in this work we define k to be in the z-direction, so that 
k = (0, 0, k z ), and we only need to consider [i, v € {x, y}. 



From symmetry arguments^ we know that C"" is non-zero 



only if fi = v, and C * is non-zero only when fj, ^ v. Also, 
= when 9 = ttM 
In our Monte Carlo we have access to the variables 



J 2 and 



In order to monitor all correlators involv- 



ing the Ji variables, we need to write C^(fc) and C}£(k) 
in terms of the correlators that we can measure: C^Ak), 
(QinifyQi^—k)) and {Ji^fyQi^-k)) . It is easy to ar- 
gue that (J 2A1 (fc)<5i Al (-fc)) = (J 2 , l (-fc)Qi /i (fc)), which are 
the only non-zero cross-correlators of J 2 and Q%. To obtain 
expressions for C^(fc) and Cj£(k) we can couple the orig- 
inal J variables to external probe fields A ext by adding the 
following terms to Eq. ((TJ: 



5S = iJ2 [M-k) ■ A xt (k) + J 2 (-k) ■ A c 2 xt (k) 



(13) 



We carry the fields through the duality procedure which 
leads to Eq. (|8). By taking derivatives of the resulting partition 
sum with respect to the external fields, we can derive expres- 
sions for Czz and C 12 v in terms of correlators which we can 
measure: 



^ll(k) — 



1 (27T) 2 (Q ly (k)Q ly (-k)) 
Mk) \.h\ 2 vi(k) 2 



(14) 



9 2 C 22 y (k) 2(2n)0(Q ly (k)J 2 y(-k)) 



|/k| 3 «i(*0 2 |/fe| 2 «i(fc) 2 

[2n(Q lv (k)J 2v (-k))+eC^(k)] 
2 sin ■ vi (k) 



Cll{k) = 



(15) 



To be explicit, in the above equations we have set [i = x, 
v = y. We note that on the symmetric line v±(k) = v 2 (k) 
and so C^(fc) = C^ 2 (k). Whenever we present numerical 
data on the symmetric line, we have performed appropriate 
averages over both of these measurements and all directions 
to improve statistics. 

In order to determine the critical exponents of the model 
at various phase transitions, we will also monitor the deriva- 
tives of the correlation functions with respect to parameters 
in the potential. One option is to study derivatives with re- 
spect to g (here we are working on the symmetric line where 
9i = g 2 = <?)■ However since g controls marginally long- 
ranged interactions in our model, it is possible that universal 
properties, and in particular critical exponents, might depend 
on iti^i To avoid possible difficulties in interpretation due to 
driving the transition while varying g, we have chosen to in- 
troduce a short-range interaction into the potential, so that 



Mk) = ^+ ti, 

\jk\ 



U2(fc) = 77T + t2 ' 

\Jk\ 



(16) 



where t\ and t 2 are parameters controlling the strength of the 
additional short-range interaction. We can drive transitions by 
varying t\ and i 2 , with the expectation that critical indices will 
depend only on g and ?/. We can fix g at its critical value g cr it, 
which we will find using our modular group analysis. We will 
extract critical exponents by taking derivatives with respect to 
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t\ and t%, at t\ — <2 = and g = g„it. We will see in Sec. I VII 
that we need symmetric and antisymmetric combinations of 
f i and £2 to extract the critical exponents. We define t s to be 
the short-ranged parameter in the symmetric direction, and t a 
in the antisymmetric direction, which leads to t 1 / 2 =t s ±t a . 
When computing the derivative of a general expectation value 
(O), of an observable O, we use the following formula: 



d(0) 



dt 



s/a 



(O) 



,=0 



dS 



dt 



s/a 



o 



dS 



dt 



sj a 



(17) 



The action S is the action given in Eq. ((H), which is what is 
used in the Monte Carlo. 

The current-current correlations C"* represent the response 
of the current J a ^ to an externally applied field Af^. We can 
view our system with long-range interactions as having an- 
other, internal, gauge field, induced by the other currents in 
the system i 21,23 - 24 In systems with short-range interactions, the 
quantity C^(fc) • L, with k = k m - in = (0, 0, can be used 
to detect the phases of the system because it decreases with 
system size L when the J variables are gapped and increases 
when the J variables are condensed. This allows the loca- 
tion of phase transitions to be determined by finding crossings 
°f Cxx(^min) ■ L at different L. However, the long-range in- 
teractions in our system prevent C^(fc m m) • L from increas- 
ing when the J variables condense ^ 14 ' 21 ' 24,25 so we cannot use 
crossings in this quantity to locate the phase transitions. To 
solve this problem, we study "irreducible responses", which 
measure the response of J to the total field made up of A cxt 
and the internal field. These responses are related to the con- 
ductivities of the system. The derivation of these responses is 
given in Ref. Qj], and the result is the following equation for 
the conductivities: 



V = 



— C(l- VC) _1 , 

IM 

Wl(fc) 



cr = 



(18) 



2sin(fe 2 /2) 
2Bin7fcl/2) U 2(fc) 



CT^(fc) relates the current induced in the J\ variables in the 
x direction to the total field in the x direction, coupled to the 
same variables. crl y (k) relates the current induced in the J\ 
variables in the x direction to the total field coupled to the J2 
variables in the y direction. In Ref. we showed that con- 
ductivities such as a\ y or a^ x are zero in our system. When 
we present numerical data we take appropriate averages over 
both species and all directions to improve statistics. Unlike 
the current-current correlators, such tr^(fc m i n ) increase with 
L in the phase where the Ji and J2 variables condense, even 
in the presence of long-range interactions, and therefore this 
quantity can be used to determine the phase transitions. 



m. PHASE DIAGRAM OF THE MODULAR INVARIANT 
MODEL 

We now wish to use the modular invariance of our ac- 
tion to determine the phase diagram of the system with the 
J\ <-> J2 interchange symmetry, in the phase space defined 
by the intraspecies interaction g and the statistical interaction 
77 (r) = 9/2tt). To begin, consider the action in terms of the J 
variables given in Eq. ((T), using the potential in Eq. (f2) with 
gi = 92 = 9- The behavior of the J variables is determined by 
the parameters g and 77. We know that as g — > 00, the system 
will be in phase 00, where the J variables are gapped. As g 
decreases, the J variables will condense. To find the location 
of the phase transition, consider the action after the applica- 
tion of the duality operation S. This action is in terms of the Q 
variables. Due to the fact that V(r) ~ 1/r 2 , the Q variables 
have the same kind of interaction as the original J variables, 
with parameters gduai, ?7duah given by Eq. (|6), instead of g, 77. 

Consider the model at 77 = 0. In this model, the two species 
of loops are decoupled, and #duai = 1/5- There are two 
phases, one phase with the J variables gapped (which we call 
phase 00) and the other with the Q variables gapped (phase 
0). The phase transition between these two phases must occur 
at g = ffduai = 1- Such single loop models with long -ranged 
interactions were studied in Refs.l2llandl26[ 

Next, we can see from Eq. ((1) that our model at phase space 
coordinates (77, g) is mathematically equivalent to the model 
at (—77, g), after making the change of variables 3\ — > —J\, 
while leaving J2 unchanged. Therefore, away from 77 = 0, 
we can use the equivalence between 77 and —77 to see that 
the phase transition will again occur when g = gdu&U which 
means that near 77 = the transition between phase 00 and 
phase occurs when g 2 + rj 2 = 1. We conjecture that this 
is the case for — | < 77 < |. We will see that this conjec- 
ture leads to a phase diagram which has the same properties 
after any operation by the modular group. This phase dia- 
gram is in agreement with our numerics. Therefore we be- 
lieve that the conjecture is correct. We know that the phase 
diagram is periodic under integer shifts of 77, so in the region 
— i + n < rj < i + n we expect a phase transition out of 
phase 00 when g 2 + (77 — n) 2 = 1, Phase 00 is located in the 
region of parameter space above these phase transitions. 

We can now use the duality transform to deter- 
mine that phase is located in the region where 
(?7duai(??,5)>5duai(?7,.9)) lie in the 00 parameter region, 

i- e - Sduai + (^ dual ~ n ) 2 > l ' for some n and I 7 ? ~ n l - \ ( see 
also Fig. |2). 

When deriving the extent of phase 0, we performed the fol- 
lowing steps. First, we applied an operation of the modular 
group (specifically, duality) to the action in Eq. ((T). This gave 
us an action in terms of new variables (specifically, we ob- 
tained the Q variables). The new parameters in the action, 
Sduai and ?7duah were functions of the original parameters g 
and 77. Note that both actions refer to the system at a sin- 
gle point on the phase diagram. By determining which (77, g) 
gave (?7duai, ffduai) in the 00 parameter region, we were able 
to determine the extent of phase 0, where the Q variables were 
gapped. 
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We now want to generalize this procedure to everywhere 
in the phase diagram. This requires us to apply modular 
group operations more complicated than duality to the action 
in Eq. ([T), and so we must determine the new phase space 
coordinates (fj, g) that result from a given modular group op- 
eration. To do this we combine Eqs. and (O to get: 511 



(19) 



d 



where a,b,c,d are integers and ad — be = 1. To find the 
a, b, c, d that correspond to a given set of S and T, we write 
them in matrix form: 



a b 
c d 



and the operations can also be represented by matrices: 



S = 

fjnn 



-1 

1 

1 n 
1 



(20) 

(21) 
(22) 



We can find the a, b, c, d that correspond to a given operation 
by multiplying these matrices. Such matrices where A and 
—A are identified make up the group PSL(2, Z), which is 
equivalent to the modular group. 

We know that there is a phase transition at g 2 + t] 2 = 1, 
for 1 77 1 < |. If we apply a modular group operation, we 
will obtain an action in terms of variables with interactions 
g, fj, and these variables will have a phase transition whenever 
g 2 + ff = 1, \ fj\ < \- Therefore we can find all of the phase 
transitions in the diagram by finding all the different values 
of g, i] which have this property for some modular operation. 
The resulting phase diagram is shown in Fig. [2] We have only 
shown one period of the phase diagram, with < r\ < 1, but 
the same structure repeats for all 77. As the modular group is 
an infinite group, there are an infinite number of phase tran- 
sitions, and so our diagram does not show all of the details at 
small g. 

The solid symbols in Fig.|2]show the locations of the phase 
transitions determined numerically. In a physical system 77 
must be fixed, and so we took data in sweeps at fixed 77, vary- 
ing g. This corresponds to vertical lines in the phase diagram. 
We determined the locations of phase transitions by observ- 
ing peaks in the specific heat. We have also observed that 
a xx(^tnm) diverges with system size in phase but decreases 
with system size in all of the phases neighboring phase (this 
will be explained in Sec. lIVI ). Therefore we were also able to 
use crossings in this quantity to locate the phase transitions. - 

Finally, using Eqs. ( TT4b and ( Tl3T > we can determine that at 
the oo-O phase transition, where the J variables and Q vari- 
ables see the same potential (but opposite statistical angle), 



[gC 1 x i(k)- V Cll(k) S gn(k z )]= l M. 



(23) 



In addition, we find that in the thermodynamic limit the above 
quantity divided by \f^ \ approaches zero in phase 00 and a dif- 
ferent finite value in phase 0, and so we used crossings of this 
quantity as another way to find the location of the transition. 
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FIG. 2: The phase diagram for the "symmetric" model with <?i = 
Q2 = <?, for one period in 77. The J variables are gapped in phase 
00, and the (dual vortex) Q variables are gapped in phase 0. In other 
phases, the gapped particles are linear combinations of J and Q. Ev- 
ery phase can be mapped to phase 00 by an operation in the modular 
group. Solid symbols show where the locations of phase transitions 
have been confirmed by our numerical study. There are infinitely 
many phase transitions in the model, so at small g our diagram does 
not show every transition. The labels on the phases are explained in 
the text. 



We have mentioned above that by applying an operation 
of the modular group we can express the original problem in 
terms of new loop variables. For a given phase, if we choose 
the modular group operation which gives (fj, g) in the 00 pa- 
rameter region, these loop variables will be gapped quasipar- 
ticles in that phase. In the following we will show how to 
determine the precise physical nature of the phase and these 
quasiparticles. Starting with Eq. <£TJ, we perform a duality 
transformation on the J\ variables to obtain Eq. ([SJ. We then 
make the following substitutions: 



Gi = cQi - dJ 2 . 
F 2 = aQi - bJ 2 , 



(24) 
(25) 



A change of variables like the one above will not always map 
the independent, integer-valued variables Q±, J 2 to new inde- 
pendent, integer valued variables. If a, b, c, d are allowed to be 
any integer, G\ and F 2 may not be independent and therefore 
the action in terms of these variables will not be equivalent to 
our original action. However, if a, b, c, d represent an element 
of the modular group, then the above substitutions represent a 
valid change of variables. Note that for the duality transform 
(a, b, c, d) = (0,-1, 1, 0), and so this transformation gives us 
G\ = Q\, F 2 = J 2 . 

After performing the above change of variables, we arrive 
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at the following action: 



S[G!,F 2 ] 



E 



(2tt) 



1 



2 t IAI 



-l^2(fc)| 2 



'/ 2 



+ ^ ) |GiW| 2 -^2(-fc)-Gi(fe) 



where 



.9 



(d + tic) 2 + g 2 c 2 ' 
(b + na) (d + rjc) + g 2 ca 
(d + rjc) 2 + g 2 c 2 



(26) 

(27) 
(28) 



In the above we have specialized to the potential given in 
Eq. (O, on the symmetry line where i>i = v 2 - One can 
check that the new parameters fj, g are precisely those given 
by Eq. (fT9l l, We can then dualize the F 2 to obtain the follow- 
ing action: 



S[Gi , G 2 



1 2ng 

2 r \fk\ 



^z27T7?Gi(-fc) • aG 2 (k), 



(29) 



with G 2 = V x ac 2 ■ We can now understand the gapped vari- 
ables as linear combinations of J variables and Q variables. 
A more general version of the above equation is given in the 
Appendix. 

In the above we chose a modular operation to map a given 



77, g to the region of the phase diagram where g 



+ Tf > 



1) \ fj\ — \- We could instead have chosen a different oper- 
ation which mapped to g 2 + (fj — n) 2 > 1, \ fj — n\ < h, since 
this would still be in the 00 parameter region. The coefficients 
of this new transformation would change by a —> a + nc, 
b —> b + nd, with c and d remaining unchanged. We can see 
from Eq. ( T2~5l ) that the new variables G\ (and by symmetry 
G2) are the same regardless of which part of region 00 the 
original model is being mapped to. In what follows we will 
always choose the modular operation which maps to the re- 
gion g 2 + fj 2 > l,\fj\ < i. We have found that all physical 
results depend only on the coefficients c and d, which are the 
same regardless of which part of the 00 parameter region the 
g, fj variables are in. 

We will label each phase by the label — . (Note that for 
phases with < 77 < 1, c and d have opposite signs). For any 
modular transformation, c and d are mutually prime, so this 
label will be an irreducible fraction which uniquely identifies 
c and d. This label is practical for a number of reasons. From 
Eq. (l25T l we see that it gives the nature of the gapped quasi- 
particles in this phase. It also gives the 77 value at which this 
phase touches the g = axis, which is also the 77 value which 
maps to g = 00. Phase 00 with c = 0, d = 1, and phase 
with c = 1, d = both conform to this label. Figure[3]shows a 
section of the modular invariant phase diagram with the labels 
assigned. 

We can understand each phase as a condensate of objects 
which have G\ = or G 2 = 0. An example of such an 



object would have Q\ — d, J 2 = c. In our Monte Carlo 
simulations we can greatly reduce autocorrelation times by 
attempting worms of these composite objects. 

Let us provide some examples of the application of the 
above approach. Consider performing an experiment on this 
system by decreasing g while holding 77 constant at 7/ = -, 
with n an integer, and n 7^ 2. This is equivalent to a vertical 
sweep in Fig. [2] or a sweep along the symmetric line in a fig- 
ure similar to Fig. Q] At large g, the J variables are gapped. 
As g is decreased the J variables condense and the Q vari- 
ables are gapped. In fact, the precise meaning of condensation 
of the J variables is that their dual Q variables are gapped^ 
Though the intraspecies potential of the J particles (which is 
controlled by g) is small, the J variables feel a statistical in- 
teraction, and this limits how many loops of the J variables 
can form. Also, g is large enough that large values of J are 
still costly, and so we expect that when the J variables first 
condense they have strength one. This means that complex 
composite objects of the J variables are not condensed here. 
Consider further decreasing g. In Eq. ©, we see that for non- 
zero 77, at small enough g the parameter gduai will also become 
small. Now we enter a phase where both the Q and J variables 
want to proliferate in some form. The gapped variables are 
now the G variables given by the modular operation ST n S, 
which has coefficients (a, b, c,d) = (—1, 0, n, —1), leading to 



1 

n 2 g 



n 



1 

n 



(30) 



As g — > 0, these variables will have g — > 00, and therefore 
no loops. The small value of g thus leads to the binding and 
condensation of more complicated composites of J (in partic- 
ular Qi = 1, J 2 = —n). These objects will see no statistical 
interaction, and loops of these variables can form more easily 
than loops of the J variables in phase 0. Specifically, under 
the change of variables in Eqs. ( f24b and ( |25| ). the interactions 
in Eq. d26l i are such that the G\ variables want to be gapped 
and the F 2 variables condensed (hence the G 2 variables are 
also gapped). This is phase "1/n". The transition from phase 
to phase 1/n occurs at g = \jn(\/n 2 — 1). 

Now consider the same experiment as above, this time hold- 
ing 77 = |, as shown by the vertical line in Fig. [3] Phase 
00 and phase have the same properties as in the earlier 
case. At g = g, the system enters phase 1/2. The new 
gapped variables in this phase are related to the J variables 
by the operation ST 2 S, and so have Gi = 2Qi + J 2 - At 
g = 1/ (5\/21) the system enters phase 2/5. The new gapped 
variables correspond to the operation ST~ 2 ST 2 S, which has 
(a,b,c,d) = (—2,1,-5,2). They remain gapped even as 
g — > 0. The new condensed variables see no statistical in- 
teraction and can condense completely. 

In the general case of rational 77 such that 77 = -, with r and 
s mutually prime integers, we can find a modular transforma- 
tion such that c = s, d = —r. For such a transformation, 
we can see from Eq. ( |27] > that g = l/(gc 2 ) — > 00 as g — > 0. 
Therefore for general 77 the system will pass through a num- 
ber of phases with different gapped variables, before finally 
reaching a phase on the 5 = axis where c and d are related 
to 77, and which can be viewed as a condensate of composite 
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FIG. 3: A section of Fig. [2] blown up to show the labels on the vari- 
ous phases. Not all transitions are shown as there are infinitely many 
of them at small g. The dashed line is at rj = 2/5, where the data in 
Figs. [4][5] and|6]were taken. 



objects like Qi = r, J2 = s. 



IV. CORRELATION FUNCTIONS AND CONDUCTIVITIES 

In our Monte Carlo, we can measure the correlators be- 
tween the current variables, C} 1 (fc) and C 12 (k), where the 
J subscript refers to the fact that these are correlators in the 
J variables. Here and below we are dropping the direction 
subscripts on these variables: C 11 means C]} x and C 12 means 
Cl 2 . We would like to determine the values of these correla- 
tors in the thermodynamic limit for all the phases in the phase 
diagram. The conductivities cr^(fc), <r^ 2 (k) are functions of 
these correlators, so this will also give us the values of these 
conductivities. We know the values of the correlators in phase 
00, because in this phase the J variables are gapped. This 
means that the only excitations are small loops, which im- 
plies that Cj (fc) ~ fc 2 . Since we measured the correlators at 



1 

L- 



k = (0, 0, k z = ?f) = fc min , we find that C n (fc rom ) 
The smallest excitation that contributes to C 12 (k) consists of 
a small loop in each of the J\ and J2 variables. An estimate of 
such contributions with cross-linking between the loops leads 
to Cj 2 (fcmin) ~ — sin(0)fe roin ~ 1/L 3 . From these corre- 
lators we can determine that the conductivities vanish in this 
phase. 

In the previous section, we have interpreted each phase by 
going to the appropriate G\ and G2 variables, and since these 
variables are gapped in this phase, we know the behavior of 
the G correlators for the same reasons given above. Therefore 
we wish to express C} 1 (fc) and C] 2 (fc) in terms of C^. 1 (fc) and 
C(?(k), where the latter are correlators of the new variables. 
To do this, consider the combination 



Dj{k) 



Cfjk) , iGf(k) 



siri if 



I sin % I 



(31) 



Consider the effect of the duality operator S on this object. We 
can derive the following relation between the complex correla- 
tion Dj(k) in the direct variables, and the complex correlation 
Dg(k) in the dual variables*^ 



D Q (k) = z 2 Dj{k) + z. 



(32) 



Note that this equation is a relation between two different cor- 
relation functions at the same point in the phase diagram. The 
periodicity operator T does not change the correlation func- 
tions. Combining the actions of the two operators leads to the 
following relation between correlation functions Dq (fc) of the 
G variables in Eq. d3Tl i. and the correlation functions Dj(k) 
of the original J variables;^ 



D G (k) = (cz + d) 2 Dj{k) + c(cz + d), 



(33) 



where c and d are the parameters of the modular group oper- 
ation which gives the gapped quasiparticles. We can rewrite 
Eq. d33l to get expressions for the J correlation functions in 
terms of only the G variables: 



C)\k) 





eg 



[{eg + d) 2 - c 2 g 2 ]C}}{k) - 2cg(cr 1 + d) sgn (k z )C^?{k) 



7T (crj + d) 2 + c 2 g' 2 



[(cr) + d) 2 + c 2 g 



212 



(34) 



, sin 
Cf(k) 



-c(cr] + d) [(cr) + df - c 2 g 2 ]C 12 (k) + 2cg(cr 1 + d) sgn (fc z )C n (fc) 



7T (cry + d) 2 + c 2 " 2 



[(ci] + d) 2 + c 2 g 2 } 2 



(35) 



In the phase where G are gapped, Cg^fcmin) ~ 1/L 2 and 
Cq 2 (fc m in) ~ 1/L 3 , and so in the thermodynamic limit the 
behavior of C} 1 (fc m i n ) and C} 2 (fc m i n ) is given by the first 
terms in the above expressions. Our numerical results agree 
with this analysis. Plots of these first terms compared to the 



numerical data are given in Fig.|4]for Cy-(k m i n ), and Fig. [5] 
for Cj 2 (fc m i n ). To find the curves that correspond to the the- 
oretical predictions, one reads off c and d from the label of a 
given phase, and substitutes these coefficients into the leading 
terms of Eqs. (f34t and d35l l. 
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FIG. 4: C} 1 (fcmin) ■ L as a function of g for r\ — | (i.e. along the 
dashed line in Fig. [3), for system sizes L — 6 and L = 8. The 
dashed curves correspond to the theoretical predictions in the differ- 
ent phases. The dotted vertical lines denote the phase transitions, 
which occur at g = | and g = - \— . 



From the correlation functions we can also determine the 
conductivities. We find that in the thermodynamic limit (and 
assuming d ^ 0) 



°xa:(fcmin) 0, 
^ (fcmin) - 



(36) 
(37) 



These conductivities are determined solely by the coefficients 
c and d, and hence by the "— d/c" label. Figure [6] shows 
<T%y(k m in) for V = |. As g is decreased, the system passes 
through both phase 1/2 and 2/5, and the conductivity takes the 
expected values of — = ^2 and ^| in these phases. At 
a phase transition the G variables are not gapped and so the 
above expressions do not hold. In phase 0, d = 0, and so to de- 
termine the behavior of the conductivities one must take into 
account the subleading terms in Eqs. (T3~4l and ( 1351 1. When this 
is done one finds that CT^(fc m j n ) diverges in phase (which is 
why its crossings can be used to detect phase transitions) and 
a^y (fcmin) approaches a non-universal value. 

Note that in the above expressions, if we chose a differ- 
ent modular operation which had fj — >• fj + n, this would not 
change Dq(U), nor would it change the modular group co- 
efficients c and d, and therefore the above equations would 
remain unchanged. Therefore shifting f] by an integer, which 
is equivalent to choosing a different modular operation to de- 
scribe the gapped particles, does not change any of the physi- 
cal properties. 

A different situation arises when we discuss shifting ij by 
an integer, i.e. 77 — 77 + n. This corresponds to describing 
a different point on the phase diagram, for example a point 
with i) » 1/3 after a shift of 1 would have 77 w 4/3. From 
Eq. <Q}, it is clear that the correlators Cj l {k), C) 2 {k) should 
have the same properties after the shift, but this is not obvious 
from Eqs. (|3~4l and ( f33T >. However, in order to get an action in 
terms of gapped quasiparticles in a phase at the shifted rj, we 
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FIG. 5: Cj (fc m i n ) • L as a function of g for the same system as in 
Fig. E] The dashed curves correspond to the theoretical predictions. 
The dotted vertical lines denote the phase transitions. 



T- X 

b 

CM 




FIG. 6: 2iraly(k m i n ) as a function of g for the same system as in 
Fig. [4] We observe that the conductivity approaches ^2 in phase 
1/2 and in phase 2/5, as expected. The dotted vertical lines de- 
note the phase transitions. In phase 0, a^y approaches non-universal 
values. 



must apply a T _n operation to our action before we apply the 
modular operation for unshifted rj. This changes the modular 
coefficients b and d: b —> b — an, d — s- d— cn, and this cancels 
the shift in 77 in Eqs. ( t34l > and ( t35T > to leave the correlators un- 
changed. Therefore when r\ is shifted different quasiparticles 
become gapped. 

Though the current-current correlators do not change when 
T] is shifted, the conductivities do change [note the dependence 
on d in Eq. (l37lil. Though the equivalence of the correlators 
implies that the system's response to an applied field is un- 
changed by a shift in 77, this does not mean that the system's 
response to the internal fields is unchanged. In particular^ 
when defining the conductivities we are treating the statisti- 
cal interaction as a long-ranged interaction mediated by real- 
valued internal gauge fields. Shifting 77 changes the strength 
of this interaction, which in turn changes the action of the in- 
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ternal gauge fields. This is responsible for the change in the 
conductivity. An interesting case is the effect of such a shift 
on the conductivity in phase 0. In this phase the J variables 
are condensed, and since these are the variables which carry 
the current in this phase cr** diverges, as we have seen. We 
can apply the operator T 1 to phase to get phase 1 . In this 
phase, the partition function for the J variables is exactly the 
same, but cr** does not diverge. To understand this, recall the 
precise meaning of condensation: a variable is said to be con- 
densed if its dual variables under the formal duality transfor- 
mation are gapped. The variables dual to the J variables are 
the Q variables that are gapped in phase 0. However, the Q 
variables are not gapped in phase 1 and hence the J variables 
are not condensed in the above sense. Instead, some other 
variables, which can be determined from the substitutions in 
Eqs. (l24t and d25l l appropriate for phase 1, are condensed. An- 
other way of interpreting condensation is that in calculations 
like the current-current correlations we can replace integer- 
valued condensed variables by real-valued variables, and per- 
form Gaussian integrals over these variables. By the above ar- 
guments we can do this in phase but not phase 1, and noting 
how we defined the conductivity for the J currents, it "knows" 
whether or not this real-valued replacement is possible. This 
explains the difference in conductivities between phase and 
phase 1. 

The phase diagram in Fig.|2]has a number of special "fixed" 
points which are unchanged by an operation of an element 
of the modular group. There are two types of such points: 
"triple points" where three phase transitions meet, and points 
halfway along a phase transition line, such as the point at 
i] = 0,g = li^ii The invariance under a modular operation 
means that in Eqs. d34]i and (O we have C} 1 ^) = Cg(fc) 
and C} 2 (fc) = C^?(k). We can then solve the two equations 
to determine the correlation functions, and therefore also the 
conductivities. We obtain the following results, applicable at 
all fixed points with g > 0: 



Cf{k) 



sin(4f)| 
2irg 



cy(k) = o. 



We can then determine the conductivities: 



2ir{g 2 + t] 2 )'' 
V 

2Tt(g 2 + ^Y 



(38) 
(39) 

(40) 
(41) 



We have verified these equations numerically for the fol- 
lowing points (rj, g): (1/2, ^3/2), (1/2, 1/2), (1/2, y/3/6). 
When i] = the two species of particles are decoupled, and 
we studied this system in Ref. |2lJ, We found the above equa- 
tions to hold for the fixed point (r/ = 0, g = 1). 



V. NATURE OF TRANSITIONS 

Our numerical approach allows us to study the properties 
of the various phase transitions in Fig. [2] We have attempted 




FIG. 7: Histograms of e at various points on the phase transition 
between phase oo and phase 0, using sizes L — 10, 12, 14, 16. A 
normally-distributed histogram implies that the transition is second- 
order, while a transition with multiple peaks implies first-order. The 
first panel corresponds to r\ = \, the next two panels correspond to 



V 



We see no evidence of first-order behavior at these sizes. 



to determine the order of the transition between phase oo and 
phase 0. To do this we study histograms of the total energy 
e at the phase transition, utilizing the fact that we know the 
exact location of the transition. In a second-order transition, 
we would expect such histograms to be normally distributed, 
while for a first-order transition we may see multiple peaks 
in the distribution. Figure [7] shows histograms taken on the 
oo-O transition, at 7/ = 4, j, §. Histograms for system sizes 
L = 10, 12, 14 and 16 are shown. We see normally dis- 
tributed histograms, suggesting a second-order transition. We 
can show that in our sign-free reformulation using the Qi, J2 
variables, (e) = 1 — for the model with J\ -H* J2 inter- 
change symmetry, at all values of g and 77. Our Monte Carlo 
measurements of (e) confirm this. 

The modular invariance of the system implies that all phase 
transitions that are related by a modular operation will have 
the same properties. In fact, one can show that in our vari- 
ables Qi, J2, any two points related by the modular group 
produce simulations with the same energies, so the histograms 
should be identical. In these variables the updates used in 
the Monte Carlo are different, but if they are done properly 
the results should be the same. We will check this by study- 
ing the properties of the line of phase transitions at 77 = i. 
There are two modular group operations which map the 00- 
phase transition to this one. The first is T 1 ST 1 S', which 
has (a, b, c, d) = (0,-1, 1, —1). This maps the three above 
points at 77 = g j 3 j 5 to three points with 77 = | and 
g = s/Q/4, \/l5/6, 1/vZ Histograms at these points are 
shown in Fig. [8] Once again, we see no evidence of first- 
order behavior at these system sizes. The second modular 
group operation which maps the oo-O transition to the line at 
?7 = 2 is ST- 1 S, which has (a, b, c, d) = (-1, 0, -1, -1). 
This maps the three points on the semi-circle to 7/ = i and 
g = 1/ v / 6,a/15/10,1/(2V5). Histo grams for these points 
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FIG. 8: Histograms of e at various points on line of phase transitions 
at rj = |, using sizes L = 10, 12, 14, 16. Each point is related to 
a point on the oo-O phase transition by the operation T 1 5T S. The 
first panel maps to g = g, the next two panels map to g = 
We see no evidence of first-order behavior. The histograms are also 
identical to those in Fig. [7] which provides a check on our Monte 
Carlo. 
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FIG. 9: Same as Fig. [8] but the operation to map these points to oo-O 
phase transition is ST~ 1 S, as explained in the text. 



are shown in Fig. [9] and they also show no sign of first-order 
behavior. The histograms for the related points in Figs. |7l [S] 
and [9] are identical, as predicted by the above argument. 

We have also studied the system at the "triple points" on 
the modular invariant phase diagram, where three phase tran- 
sitions meet. We expect all such points to have the same 
properties, and we have studied the points at the ends of the 
7] = h line of phase transitions, which occur at g = \[2>j2 
and g = \/3/6. Histograms for these points are shown in 
Fig. [Tol a) and (b). We see that the histograms have two clear 
peaks, indicating that these are first-order transitions. 

In Ref . I2TJ we studied the phase transition at (77 = 0, g = 1) 
and found it to be continuous. This point maps to the point 
77 = i, g = -i, and we have studied the phase transition at this 
point to confirm the second-order behavior. Histograms at this 
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FIG. 10: Histograms of e at some special "fixed" points in Fig. [2] 
using sizes L = 10, 12, 14, 16. The first two panels show histograms 
at "triple points" where three phase transitions meet. Panel (a) shows 
the triple point at the upper end of the line of phase transitions at 
g — i, which has g — \/H/2. Panel (b) shows the point at the 
lower end with g = V3/6. Both histograms have a double-peaked 
structure which indicates that the transitions at these points are first- 
order. Panel (c) shows histograms at the fixed point g = \,g = 
I . These histograms show no sign of first-order behavior, which we 
expect since this point should have the same properties as the point 
g = 0, g = 1 which is known to be continuous.— 

point are shown in Fig.fTol'c). and we see no sign of first-order 
behavior. 



VI. CRITICAL EXPONENTS OF PUTATIVE 
SECOND-ORDER TRANSITIONS 

Apart from the triple points, the transitions we have stud- 
ied are second-order, and we can now determine their critical 
exponents. In Fig. [2] the oo-O transition seems to take place 
at an ordinary critical point, however in Fig.[JJwe see that the 
transition actually takes place at a tetracritical point when we 
allow V\ and v-i to be not equal. Figure [TTT a) shows a closer 
look at the upper tetracritical point in Fig.[JJ At such a tetra- 
critical point there are two scaling directions, each with a dif- 
ferent correlation length exponent^ Our system is symmetric 
under the interchange of J\ and J 2, which implies that one 
scaling direction is in the symmetric (8v\ = 61)2) direction, 
with critical exponent v a , and the other is in the antisymmet- 
ric direction (8v\ = —8V2), with critical exponent v a . We can 
determine which of these exponents is larger based on how the 
phase transition lines meet. A simple renormalization group 
argument shows that the phase boundaries in the local coordi- 
nates A s , A a obey 

A ~A?/"-. (42) 

Our data shows that the phase transition lines have the same 
shape as those in Fig. QT| and this combined with the above 
equation implies v s > v a . 
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FIG. 11: (a) Schematic blow-up of the phase diagram in the V\ , V2 
(Fig. Q} variables, near a tetracritical point. At such a point we ex- 
pect two scaling directions with distinct critical exponents v a and v a . 
Due to the symmetry on interchange of Ji and J2, we expect the 
scaling directions to be in the symmetric and antisymmetric direc- 
tions, shown by the dotted lines. The shape of the phase transition 
lines implies that v s > v a . (b) Phase diagram for the special case of 
6 = tv (77 = 1/2). We can see that along the symmetric direction 
the system goes along the phase boundary, and we can only drive a 
phase transition across the antisymmetric direction. 



We will extract the critical exponents by taking appropriate 
derivatives of C 11 (fc m ; n ) ■ L. In this system with long-ranged 
interactions, C 11 (fc m i n ) ■ L approaches a constant (possibly 
zero) value in each phase. At a critical point, it jumps from 
one value to another, which leads to a peak in its derivative. As 
mentioned in Sec. [II] we want to take derivatives with respect 
to a short-ranged part of the potentials, given by the param- 
eters ti and t2 in Eq. dl61 l. To extract v s we take derivatives 
with respect to the symmetric combination t s , while for v a we 
use the antisymmetric combination t a . C 11 (fc m i n ) • L has the 
scaling form: 

C 11 (fc m ; n ) • L = f a (Lt" 3 ) symmetric direction, 
C 11 (fc m ; n ) • L = f s {Lt v a a ) antisymmetric direction, 

where f a and f s are scaling functions. This leads to 



^^J 1 (^min) ' L 



dt s 

dcy(k min )-L 



dt a 



-l/Vs 



(43) 
(44) 



so v SXL can be extracted by fitting curves of d(CY(k m i n ) ■ 
L)/dt s a vs L. Such curves, at the oo-O transition, are shown 
in Fig. \l2\ a) for the symmetric derivative and Fig. \l7\ b) for 
the antisymmetric derivative. The extracted values of v s and 
v a are shown in Fig. Q~3] We see that v s > v a , as expected 
from the shapes of the phase transition lines near the tetracrit- 
ical point. We also see that neither exponent is close to 1/3, 
supporting the conclusion that we have second-order transi- 
tions at these points. The exponent v a is decreasing as we 
move along the phase transition away from 77 = 0. On the 
other hand, the error bars on v s are too large to tell whether it 
is varying. Error bars for v can be estimated from the fits to 



the d(Cy (fc m in) • L)/dt curves. However, we may have sig- 
nificant finite size effects in our results, even though we are 
simulating exactly at the transition, because we do not know 
the subleading corrections to Eqs. (l43l and f44l . To account 
for this in the error bars we performed fits both including and 
not including the data at L = 6, and if the errors from the fit- 
ting procedure were not large enough to encompass both val- 
ues we increased the error bars. The values of v were taken 
from the fits which included the L = 6 data. We have also 
plotted the critical exponents measured in Ref. 21 at the point 
g = 1,1] = 0. At this point the two species of particles are 
decoupled, so v s = v a . 

We could not determine v s on the 77 = i line because 
changing t s does not move the system through a phase tran- 
sition, instead it moves along the line of phase transitions, as 
seen in Fig. [TU b). On the other hand, we can argue that the 
transition driven by t a is equivalent to that on the oo-O line 
at the related point, and the v a values on this transition are 
shown in Fig.fOl'b). 

The transition at 77 = 5 is a transition between phase I, 
where the J\ variables are condensed and the J2 are gapped, 
and phase II where J\ is gapped and J2 condensed. The ir- 
statistical interactions prevent both types of loops from con- 
densing simultaneously. The two species could behave as im- 
miscible fluids and phase separate, or they could coexist in a 
critical soup.— ^ Our result that the transition is second-order 
implies that the two species can indeed form such a critical 
state. 

An open question is how three transitions meet at a triple 
point in the modular phase diagram in Fig. [2] All three tran- 
sitions could be second-order all the way to the triple point, 
or they could have bicritical points where they become first- 
order. Our results show that on at least part of the phase 
boundaries the transition is second-order. A more detailed 
study could determine whether the phase transitions do be- 
come first-order at some point, and where this point is. 



VII. DISCUSSION 

We studied a model of two species of particles with mu- 
tual statistics and long-ranged interactions such that the model 
has the same form after the application of operations from the 
modular group. Using this modular invariance, we were able 
to analytically conjecture the phase diagram and determine 
the values of the current-current correlations and conductiv- 
ities in each phase and at all points that are invariant under 
the modular group. We can also describe each phase in terms 
of particles gapped in that phase. Using a reformulation of 
the model that does not have the sign problem, we performed 
Monte Carlo studies and firmly established the conjectured 
phase diagram. Furthermore, we numerically determined the 
order of the transitions in the phase diagram of this modular 
invariant system. We found the triple points to be first-order 
but the other phase transitions we studied were second order. 
The second-order transitions are evidence for a novel critical 
loop soup state in the case 9 = tt. 

Exact results derived from the modular invariance greatly 
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FIG. 12: Plots of the derivatives 
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a ) ■ L)/dt 3 and 

^(Cj 1 (fcmin) • i) / 3t a at various points on the oo-O phase transition. 
Error bars were obtained by comparing the results of independent 
simulations. We expect such plots to scale as L 1 ^ 3 -". The values 
shown in Fig. Q~3] were extracted by fitting these curves to the func- 
tion a + bL 1/u . 



improved our numerical study. We were able to derive a num- 
ber of useful checks on the Monte Carlo, as well as run sim- 
ulations exactly at the location of the critical points, greatly 
simplifying our measurements of critical exponents. 

Since we have only studied this model at relatively small 
system sizes, it is possible that the transitions which we 
observe to be continuous are actually first-order. Studying 
larger system sizes could confirm the second-order behavior. 
It would also be interesting to study field theories for such 
systems with marginally long-ranged interactions and mutual 
statistics 
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Studying this model at larger system sizes 
would also allow one to obtain better estimates of v, since any 
subleading terms would have a reduced effect. In addition, at 
larger sizes one could determine the behavior near the triple 
points, in particular where the transition changes from first- 
order to second-order. However, the long-ranged nature of the 
interactions in the model means that in our reformulation, the 
amount of computer time needed scales as L 6 , making such 
larger studies difficult, but possible with future resources. 

At 9 = 7T, our (2+l)-dimensional model is relevant to the 
study of unusual phase transitions in (2+1) dimensions.^— It 
also applies to the study of (2+l)-dimensional symmetry pro- 
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tected topological (SPT) phases and phase transitions 
well as the surface states of (3+l)-dimensional SPT phases 
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FIG. 13: Critical exponents along the (a) symmetric and (b) antisym- 
metric scaling direction, extracted from the data in Fig.[T2] Error bars 
come from the fits, and are further increased to account for finite-size 
effects as discussed in the text. The blue circles represent the value of 
v obtained in Ref.[2j]. v 3 cannot be measured when 6 = tt, since the 
system does not cross a phase transition in this direction [see Fig. II II 
(b)]. On the 6 — tt line there are two points which map to each point 
on the oo-O semicircle. 



It may also be possible to use similar lattice models to study 
such SPT phases in the bulki2£d2. More generally, our loop 
model provides a precise lattice realization of a topological 
field theory. It would be interesting to study such lattice mod- 
els for other topological field theories4^^ 
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Appendix A: Modular transformations for general potentials 

The method in Sec.[nI]allows us to apply duality and peri- 
odicity operations to an action to obtain a new action in terms 
of new variables. Such a procedure is possible for any choice 
of potentials, though potentials other than the ones used in this 
work will not be modular invariant. We can still interpret the 
new variables as being gapped in a certain phase, but without 
modular invariance we cannot use this to determine the ex- 
act locations of the phase transitions, or to predict where the 
new variables will be gapped. In Ref. [HI we used numerical 
methods to determine which variables were gapped in each 
phase, and we were then able to use this procedure to find the 
action for these gapped variables. We were also able to deter- 
mine the correlation functions and conductivities in the phases 
where we knew the gapped variables, using the same methods 
as in Sec. El 

We now provide the equations that generalize the methods 
we have used in this paper to any potential. These equations 
also represent the procedure used in Ref. generalized to 
any operation of the modular group. The action in terms of 
new variables is obtained by starting with Eq. (O, making the 
substitutions Q\, J2 — > G\,F2 given in Eqs. d24b and (f25j), 



J 



and then dualizing the F2 variables to obtain the G2 variables. 
We find 



S = 



\Y.[vGAk)\Gi{k)\ 2 +v G2 (k)\G 2 (k)\ 2 

k 

^iOaGxi-k) ■ a G2 (k), 



(Al) 



where 



VG 1/2 (k) 



{2n) 2 v 2/1 (k) 



{2nd + 9c) 2 +v 1 {k)v2{k)\f k \ 2 c 2 



and 



2tt- 



{2nb + 6a) {2nd + 6c) + wi(fc)w 2 (fc)|/fc| 2 ca 
{2nd + 6c) 2 + v 1 {k)v 2 {k)\f k \ 2 c 2 



We can also express the current-current correlators in terms of 
correlators in the new variables, using the same methods as in 
Ref. [TH For simplicity, we specialize to the symmetric line 
where v\ = V2 = v: 



cY(k) 



Cf{k) 



v{k)\h\ 2 c 2 



{6c + 2nd) 2 + \f k \ 2 v{k) 2 c 2 



3c + 2nd) 2 - \f k \ 2 v{k) 2 c 2 ]C 1 G 1 (k)- 4sm^v(k)c(6c + 2nd)C}?(k) 

l(8c + 2nd) 2 + \f k \ 2 v(k) 2 c 2 } 2 {2ll) ' 



2sm!fc(6c + 2nd) [{6c + 2nd) 2 ~\f k \ 2 v{k) 2 c 2 ]C^ 2 {k)+4sm^v{k)c{6c + 2nd)C^{k) 

(2tt) . 



(0c + 2nd) 2 + \f k \ 2 v{k) 2 c 2 
When v = tj4, these equations reduce to Eqs. ( [34-b and (|35l l 



lc + 2nd) 2 + \f k \ 2 v{k) 2 c 
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